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We analyze the dynamics of structural transitions between normal right-handed B form and 
unusual left-handed Z form for a linear DNA molecule. The dynamics under the external torque 
in physiological buffer is modeled by a Langevin equation, with the potential term given by the 
authors previously [Phys. Rev. E 84, 021926 (2011)]. With this model, we first simulate the 
relaxation processes around B-form structure after sudden changes of the external torques, where 
slow relaxation ~ t~^^^ as a function of the elapsed time t is observed. Then, the dynamics of 
structural transition from Z form to B form is computed under various external torque strength. 
For small external torques, the transition proceeds via nucleation and the growth, while for higher 
torques, Z-form structure becomes unstable, and the transition mechanism is switched to a spinodal- 
like process. These numerical results are qualitatively understood by simple phenomenological 
arguments. 

PACS numbers: 87.14.gk, 82.37.Rs, 87.15. H-, 87.15. hp 



I. INTRODUCTION 

Recent advances in experimental techniques shed light 
on the response of single double-stranded DNA molecules 
to mechanical stresses, such as twisting and stretching 
[ll-Q . Various theoretical models have been developed to 
describe the mechanical responses of the molecule. For 
example, an extended theory of the classic clastic rod 
model can describe the supercoiling and its statistical 
properties Statistical models that phenomenolog- 

ically describe various structural transitions were devel- 
oped in Ref. @. Several mesoscopic models have also 
been developed, which can describe the interaction be- 
tween DNA conformation and melting structural transi- 
tion [13, In our previous paper [l^, guided by a 
gauge principle, we have developed a mesoscopic model 
of DNA that describes the interplay between the global 
configuration and the various intrinsic structures of DNA 
base pairs, such as B and Z structures. In Ref. [l^ . 
we have also constructed an effective potential that de- 
scribes B-Z transition of GC base-pair repeats, and eluci- 
dated the statistical properties that the GC repeats tend 
to make structural transition between the usual right- 
handed B form and unusual left-handed Z form under 
small external torques, in agreement with the recent ex- 
periment of Lee et al. [IJ]- The GC repeats are fre- 
quently included in DNA sequences of promoter regions, 
and furthermore, in recent studies [l^ . the Z-form struc- 
ture is related to genetic instability. 

In this paper we are concerned with the B-Z transition 
dynamics. To construct the dynamical model, we first 
consider the Lagrangian for the linear DNA molecule by 
introducing kinetic energy, in addition to the effective po- 



tential energy given in Ref. jl2| . From this Lagrangian, 
the Euler-Lagrange equations of motions are derived. To 
model the effect of surrounding physiological buffer, we 
moreover include viscous torques and random forces into 
the equations of motions. The resulting equations of mo- 
tion are utilized to analyze both the relaxation dynam- 
ics and the B-Z transition dynamics under various exter- 
nal torques. An alternative approach to B-Z transition 
was proposed in Ref. 15[ , where the structural transition 
dynamics for an idealized dissipationless DNA molecule 
were mediated by uniformly moving kinks. (See the Ap- 
pendix for the application of the alternative treatment 
to our model.) However, as we see below, the viscosity 
of the surrounding solution completely alters the transi- 
tion mechanism. In particular, the transition mechanism 
is switched from a statistically homogeneous nucleation 
process to a spatially inhomogeneous spinodal-likc dy- 
namically unstable process, as the external torque is in- 
creased. 



II. FORMULATION 

In the following, we introduce a nonlinear dynamical 
model that is a generalization of the static model given in 
Ref. To this end, we first construct the Lagrangian 
of a linear DNA molecule. The total potential energy is 
a functional of torsional angle of sugar-phosphate back- 
bones x(s) and structural order parameter p(s), where s 
is the continuous parameter representing the number of 
bases from one end of the molecule. The potential energy 
density is given by 
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where C = 91 x IQ-^oj, di = Di/tu^ with Di = 4.1 x 
10~^^ J and Wo ~ 0.6 rad/bp [l2|- The potential density 
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V{p) is given by 



V(p) = Vo[(/9M)'-l]' + W, 



(2) 



with Tc = -7.9 X 10~2iNni and Vq = 6 x lO^^Oj^ 
where the B-DNA and Z-DNA conformations correspond 
to the minima at around p ~ luq and p = — Wg, respec- 
tively. These parameter values are given in Ref. [12| to 
reproduce thermal equilibrium properties. 

The structural order parameter p(s) is related to the 
base-pair torsion, namely, the angle 6(s) subtended at 
the center axis of linear DNA, as 



dejs) 
ds 



Pis) 



Now besides the potential energy ([T|), we shall take 
account of the kinetic energy terms coming from the an- 
gular variables O and x- These are introduced in the 
following manner: The first is concerning Q. The inertia 
moment of a base pair Ii is estimated as a rigid rod of 
length L = 2nm (~ the diameter of DNA), approximately 
given by 
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1.47 X 10"*^ kg 



where Mi is the total mass of a GC base pair, 262.228 
g/mol = 4.4 X 10~^^ kg, which is estimated as the sum 
of the masses of bases G and C (151.126 and 111.102 
g/mol, respectively). Using these, the kinetic energy of 
base pairs is given by 



y \~dt 



The second is concerning x- The inertia of moment for 
sugar-phosphate backbone chain I2 is similarly given by 



AI2 X (lnm)2 = 7.8 x 10""*^ k 



g m 



where M2 is the mass of backbone chain per bp given by 
M2 = 3.9 X 10~^^ kg, which is estimated as the sum of 
two deoxyriboses (2x134.13 g/mol) and two phosphoric 
acids (2x98 g/mol). Then, the kinetic energy of sugar- 
phosphate backbone per bp is given by 



IC2 



h 
2 



dx 
dt 



In single molecule experiments, one end of a DNA 
molecule is fixed to a surface and thus, without loss of 
generosity, we impose the following condition at the end 
of s = 0: 

6(0) X(0) = 0. 

Furthermore, suppose that the external torque t is ex- 
erted on the other end s = N, where A'' is the total bp 
length of the DNA molecule. Hence, the following poten- 
tial is added to the Lagrangian: 



-TxiN). 



Putting all terms together, the total Lagrangian is given 
by 



/ (/Ci +IC2-T)ds- F„ 



L = 

= Ki + K2 +F- Fext 



(3) 



For the sake of the numerical implementation, we dis- 
cretize the argument s of dynamical variables p(s, t) and 
x(s,i), into s = 0, ds, 2ds, . . . ,nds{= N). Then, using 
Qo = 0, we have 

i 

and 

d ' 

-e,(t) = Ep,(<)ds. 



Accordingly, Ki = K-ids is given by 



K, = '-lds^± 



3 = i 



Similarly, we obtain the following discretized expressions 
for K2 and F, respectively: 
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The Euler-Lagrange equation ^ 
equation of motion for pj : 



dKi 



^ gives the 



n k 
k—i j—1 



(6) 



where fi — —dF/dpi is given by 



f, =Cds (^^-^^ Pr 

-dsV{pi). 



dids 



"^Pi - Pi-i " Pi+i 
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Similarly, from the Euler-Lagrange equation for Xii one 
obtains the following equation of motion for Xi'- 

hdsxt = 9r 

where gi = ^dF / dxi is given by 

Pi+i - Pi ( Xt+i - Xi Xi- Xi-l 



gi ~ ~ Cds 



ds 



ds^ 



ds^ 



(8) 
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Now we model the solvent effect to the DNA molecule 
by Langevin equation. For simplicity, here we neglect the 
hydrodynamic long-range interactions between DNA seg- 
ments, which would be important especially for curved, 
long DNA molecules. The Reynolds number (=inertial 
force/ viscous force) is so small in this situation that the 
viscous forces and thermal random forces thermally equi- 
librate the DNA molecule. Then, the viscous force is 
estimated as follows. When a segment of s-th bp ro- 
tates at angular velocity w = dx/dt, the solvent sticks to 
the DNA surface and the specific speed of fluid flow is 
V ~ ujR, where R is the radius of DNA. Hence, we obtain 
the velocity gradient ~ ujR/R = uj, and the correspond- 
ing shear stress ^ rjuj, where rj is the viscosity coefficient. 
In the following, the values of water (1 mPa-s) is used 
for 7]. The viscous forces applied to the DNA segment 
is given by KArjuj, where A is the area of the DNA side 
surface and a dimensionless coefficient k is introduced. 
For simplicity, A is approximated by 2Tr Rh, i.e., the area 
of side surface of cylinder with the radius R of DNA (~ 
1 nm) and the height h of the Ibp DNA segment(^ 0.34 
nm). Then, the viscous torque is given by 

(viscos force) x R ~ 2n fiR^hrjuj 
- 2.14k X 10"^°kg m^/s x uj = kuuj. (9) 

Hence, the equation of motion for Xi is modified as: 

hds^^-^ = gt - nads^ + Ci{t). (10) 
dt^ dt 

Here, the random forces ^i{t) arc introduced as Gaussian 
white processes satisfying, 

(6(0Cj(0)) = 2tiads X kBT5{t)5,j, (11) 

where fcs is the Boltzmann constant and T is the tem- 
perature of the surrounding buffer. Note that Eqs. (HU]) 
and (jlip guarantee the i-th segment's rotational energy 
l2dsx{t)'^ /2 to reach the thermal value at temperature T. 
Moreover, we postulate that the buffer does not affect the 
dynamics of p, because, in B form or Z form, base pairs 
are located inside the DNA molecule and apart from the 
buffer. 

In the limit of ds —J- 0, one obtains the following cou- 
pled stochastic partial differential equations: 

= c(M|i)-,(.,«))-VWM)) + *^. 

(12) 

and 

/a^xM) dp{s,t) \ dxis,t) , 
^^[-d7^ 

+ t6{s-N), (13) 



where ^(i ds, t) = ds is scaled random noise satisfying 

(f(si,ii)l(s2,i2)) = 2KakBTS{si - S2,ii - ^2)- 

Note here that, in the ds ^ limit with k = and 
T = 0, the resulting coupled wave equations for Q aiid y 
are formally equivalent to the model of Volkov [13, . 
These equations have a kink solution, as a special solu- 
tion, whose detailed derivation is given in the Appendix. 

III. NUMERICAL METHOD 

In this section, we develop the numerical method for 
solving ^ and ([TU| . In the following, we set ds = 1 for 
simplicity. 

First, dimensionless variables are introduced by mea- 
suring energy in units of Vo and time in units of 
^yIl/Vo{= 1.5625 X 10^^^ s). These equations read: 

n fc ,2 /,N 

EE^-^*te-x.-i-P.) 

+ dl{2pi - - pi+i) - 4[{p/uJof - l]p/wo + T*, 

(14) 

and 

(15) 

where the asterisks denote dimensionless quantities: 
C* = C/Vo, dl = di/Vo, T* = T,/Vo, T* = r/Vo, 
I2 = h/h, and a* = a/\/IiVo. Random forces £,i{t*) 
satisfy 

where the dimensionless temperature T* is introduced as 

T* = kBT/Vai^ 0.07 at 300K). (16) 

In the computation of Eq. ([T?]) . it must be solved for 
pi — d'^pi/dt*'^. Denoting the right-hand side of Eq. 
by fi, we rewrite ((T4)) as 

i n k 

Hence, 

i 

E^j' = - fi+i = ^f*- 

i=i 

We again rewrite this equation as 
j-i 
Pi + E^'j 
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and obtain the following result: 



(17) 



In our computation, Eqs. (|15p and (jl7p are numerically 
integrated. At every di-time-step computation, each ran- 
dom force in Eq. (fT5|) is set as ^idt = R, where i? is a 
random number sampled from the Gaussian distribution 
with mean zero and variance \/2Ka*T*dt. 



IV. DYNAMICS OF A LINEAR DNA 
MOLECULE 

A. Slow relaxation 

Here, we study the relaxation after a sudden change 
of external torque. To begin, we estimate two specific 
time scales: relaxation time and oscillation period. The 
specific relaxation time is given by 



relax = — = 3.6^-1 X 10-"s ^ 0.23/t-V^i/Vo 



while the time scale of oscillation is 



iosc= a/-^ = 9.258 x 10- 
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Then, at first, one might expect that the dynamics were 
classified into two: over damped dynamics if/«>l/3.6~ 
0.27; and under damped dynamics if k < 1/3.6 ^ 0.27. 
However, as we shall see below, the dynamics show more 
complicated, slow relaxations cx t~^/^, where the relax- 
ation time troiax will turn out to be the key factor for 
determining the prefactor of the slow relaxations. Here- 
after, we use dimensionless quantities and omit, for the 
sake of notational simplicity, the asterisks on the sym- 
bols. 

We simulated the relaxation dynamics at T = 0. The 
initial state was prepared as B-form structure: namely. 



UJQl 



for i = 1,2, 



where ds = 1 and the length of the DNA was set as 
n = 200. At i > 0, the external torque of t = 1 was sud- 
denly exerted on the n-th bp segment. We computed how 
the DNA state relaxed to a new stable state of twisted 
B-form structure for various dimensionless dissipation co- 
efficients K. 

Figure [ija) shows the total kinetic energies K = 
Ki + K2 as functions of time for various k. For all k 
computed, there arc slow relaxations B{K)t~^/^ , at their 
late times. The prefactor B{k) is plotted in Fig.[ljb), as a 
function of k. We clearly see the relation B{k) oc k^'^/^ 
holds. This is qualitatively understood as follows: As 
shown in the inset of Fig. [l](a), K (x t {t < 100) for the 
conserved system of k = 0. Then, the nonconservative 
systems of k ^ would depart from the behavior of the 
conserved system at the relaxation times troiax- At the 
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FIG. 1: (a) Total kinetic energies K are plotted as functions 
of time t after sudden addition of external torque r = 1, for 
K — 0.1,0.2,0.4,0.8 from top to bottom. The dashed line 
oc t"""'' is drawn for eye-guidance. Long-time dynamics of K 
for various k show slow relaxations B{K)t~^^^ . (Inset) K vs 
t for the conserved system {n — 0). (b) Plot of coefficient B 
as a function of k. The values of B are obtained by fitting 
the data of t > 10 in Fig. (a) to Bt~°-^. The dashed line is 
O.OOOOSS/t"^/^, which clearly shows B ~ k"^/^. (c) Semi-log 
plot of if vs t for the linearized system near p — uq with 



K = 0.2 and r — 1. 
exponential decay. 



The dashed line oc e 



shows the 



times, the kinetic energy K would also have their maxi- 
mum values K oc <roiax ^ K~^. After the times, K show 
slow relaxation K cx t~^/^. Putting these estimates to- 
gether, the following relation is derived: 



K 



t 



-1/2 



j-3/2^-1/2 



This qualitatively agrees with the above mentioned nu- 
merical result B{k) cx kT^I"^ . 
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FIG. 2; (Color online) Structural transition dynamics from 
Z to B form, for t — and k = 0.01. (a) Density plot 
of Pi as functions of dimensionless time t and i. (b) Semi- 
log plot of frequency distributions of B-domain sizes for t = 
799000(»), 839000(H), 879000(f) (c) Snapshot graphs of Xi 
vs i at dimensionless time t = 0, 10^, 2 x 10^, . . . , 20 x 10^ are 
plotted from below to top, respectively. 

Notice that the relaxations studied here are not usual 
exponential but an anomalously slow one ~ t^^/^. This 
type of slow relaxation is well known in the studies 
on one-dimensional nonlinear dynamical systems (e.g., 
Ref [l^ ) . For comparison, the torque response of the lin- 
earized equations near p = wo is plotted in Fig.[ljc), from 
which we see that K relaxes exponentially with t. Hence, 
the origin of the observed power-law relaxation dynamics 
is the nonlinear coupling in quasi-one-dimensional DNA 
molecules. 



B. Transition from Z form to B form 

Now we consider the dynamics of structural transition 
from Z form to B form. 

The initial state was prepared as Z-form structure: 
namely. 

Pi = -Wo, Xi = "^oi for i = 1, 2, . . . , n, 



where the length of the DNA was set as n = 200. At 
t > 0, the external torque r > was suddenly exerted on 
the n-th bp segment. We computed how the DNA state 
relaxed to the stable B-form structure from the initial 
Z-form structure, at the room temperature T = 0.07 [see 
Eq. (HID]. 

The transition dynamics from Z form to B form is 
shown in Fig. [2lja)-[2jc), for r = and k — 0.01. From 
Fig. IHa), we clearly see the following: (l)The stable B- 
form domains, indicated in light gray in the figure, are 
quickly created from the beginning of time evolution. (2) 
Then, the mixed states of B-form and Z-form structures 
are formed, where the statistical weights of B form grow 
in the course of time evolution. (3) Eventually, almost all 
segments become B-form structures. Note here that, if 
the nucleation is a random process with probability p, the 
frequency of domain size D should scale as oc . Fig- 
ure [^b) shows that this scaling relation indeed holds for 
B-domains larger than two segments and that the critical 
size of the nucleation is two segments. As a result of the 
statistically uniform nucleation, Xi changes with keeping 
the linear shape Xi ^ i [Fig. HJc)]. More preciously, the 
following approximate relation holds: 

Xi i{p) + fluctuating term, 

where (p) denotes the average value of pi. This is easily 
understood because Xi — Pj holds, and because one 
can replace pj in this equation by (p) for the statistically 
uniform nucleations. 

Now we consider how the nucleation rate depends on 
the external torque strength. When a small torque is ex- 
erted on a terminal, the interactions between neighboring 
segments will quickly induce torque balances on all the 
segments. Then, the effective potential for structural pa- 
rameter p must be modified by the exerted torque, as 
depicted in Fig. [Slja). Under a weak external torque t 
satisfying < r < T2, the activation energy for transi- 
tion from Z form to B form is given by 

AE ~ AEo - Ft, 

where AEq is the activation energy when t = 0. Accord- 
ingly, nucleation rate p must depend on r as follows: 

p cx cxp{-/3AE) cx exp(/3rr). (18) 

Figure [3Ib) plots the total numbers of B-form segments, 
Nb, as functions of time, for various r. This figure shows 
that each Nb initially increases linearly with time, which 
is consistent with the above confirmed fact that the cre- 
ations of B-form domains arc mediated via statistically 
uniform nucleation mechanism for weak external torques. 
From the fitting lines in Fig. Eljb) , we obtained the speeds 
of initial increases, dNB/dt. Figure [Sljc) plots the ini- 
tial dNB/dt as a function of t. This clearly shows the 
exponential dependency of the nucleation rate to r, as 
expressed in Eq. ^TE\i . 
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FIG. 3: (Color online) (a) Effective potentials for p are plotted 
for various external torques r = 0, ri, r2, ra, where < ri < 
T2 = 2.65 < Tz- For t = t\, the activation energy required for 
transition from Z form to B form is lowered by the amount 
indicated by the vertical arrow. Larger external torque r = T3 
makes the Z-form structure unstable. T2 is the critical value 
at zero temperature beyond which the Z-form structure is no 
longer metastable but becomes unstable, (b) Number of B- 
form segments Nb is plotted as a function of dimensionless 
time, for r = 0, 0.1, 0.2, 0.3. Thin lines indicate the linear 
fitting to the data, (c) The speeds of Nb increases are semi- 
log plotted as a function of t. 



C. Transition under high external torque 

When the external torque exceeds a critical value 
of Ti ~ 2.65, however, Z-form structure changes from 
metastable state to unstable state, as shown in Fig.|31[a). 
Hence, for such high external torque conditions, the 
structural transition would proceed via a mechanism 
other than the uniform nucleation. 

For the high external torques, the transition dynamics 
from Z-form to B-form are shown in Fig.|31[a)-|3Jc). From 
Fig. Ufa), we see that the structural transitions occur 
about 100 times faster than that for the small external 
torque cases (k = 0.01). In particular, the transition is 
not uniform: It first happens at the terminal on which the 
external torque is exerted, and then the interface between 
Z form and B form moves to the other terminal. Due to 
the B-Z interface-mediated transition, Xi makes a kink 
at the interface, as shown in Fig.SJJb). Figurc|Uc) shows 
the number of B-domain segments, as functions of 
time, where all Nb computed here are fitted well to the 
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FIG. 4: B-Z transition dynamics under high external torques, 
(a) Density plot of pi for r = 3 and k = 0.01, as func- 
tions of dimensionless time t and i. (b)xi vs i for r = 3 
and K = 0.01. The curves respectively correspond to t = 
0, 100, 200, . . . , 2000 results from below, (c) Nb are plot- 
ted as functions of t/ft for (r, k) in r = 2.5, 3, 4, 5, 6, 8 and 
K — 0.01, 0.1, 1. All data are scaled onto a master curve. 
The dashed line oc t"'^ is drawn for eye guidance (see text). 



following form 



cx 



0.5 



(19) 



This is qualitatively understood from the way of B- 
form domain growth. Suppose, at a time t, x base-pairs 
from the terminal have already made transitions into a 
B-form structure. At the moment, the Z-form segment 
adjacent to the B-Z interface is about to transition into 
the B form. Let the transition be completed within a 
time interval ttrans- Then, the sweep speed v of the B-Z 
interface is given by l/itrans- In the following, we will 
estimate ttrans- At the beginning of this transition, the 
angular velocity of the motion is assumed to be zero, be- 
cause of the large viscosity. Moreover, the B-form DNA 
of length X bp is approximated by a rigid rod of inertia 
momentum xl — x{Ii +12)- Then, its angular velocity 
r2 obeys the following equation 



xL 



dt 



XKaQ, 
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resulting in 



XKa 

Solving the following equation for itra 



2^0 



ndt. 







one obtains 



With this expression, the sweep speed is given by w = 
l/^trans- Integrating dx/dt = v ~ 1/ttrans) we obtain the 
following relation: 



4a;o Q/o kq; wqX 

T ' A 



tI 



For large x, we have 



3 tI 



tI 



t 



1/2 



and therefore Eq. (|T9l) holds. 

We have confirmed that the unstable segment adjacent 
to the B-Z interface induces the structural transition dy- 
namics. Hence, the transition dynamics can be called a 
kind of real-space spinodal decomposition, in contrast to 
the usual spionodal decomposition that occurs in recip- 
rocal space. 



V. CONCLUSION 

In this paper, we studied the dynamics of a linear DNA 
molecule. To this end, the statistical model given in 
Ref. (l2j is generalized by including rotational kinetic en- 
ergies. Then, the dynamics under the external torque 
in a surrounding physiological buffer is modeled by a 
Langevin equation. 

In Sec. lIV Al by using the Langevin dynamics, we first 
simulated the relaxation phenomena after sudden addi- 
tions of external torques. We found that the total ki- 
netic energy shows non-exponential slow relaxation as 
K = Bt-^l^, which is a characteristic of one-dimensional 
nonlinear dynamical systems. In contrast, K relaxes ex- 
ponentially with t for the linearized equations of mo- 
tions, which means that the nonlinearity is indispensable 
for the non-exponential decay. The prefactor B turned 
out to depend on viscosity proportional constant k, as 
B ~ We found that this can be qualitatively un- 

derstood from the fact that the dynamics of nonconser- 
vative systems depart from the behavior of the conserved 
system at the relaxation times of trciax ^ 1/^^. 

We then proceeded to study the dynamics of struc- 
tural transition from Z form to B form under external 



torques, r, exerted on a terminal [Sec. IIV B| . For small 
external torques, B-form domains nucleate and grow. 
This statistical process occurs uniformly within the whole 
DNA segments. The speed of B-domain growth thus 
obeys an Arrhenius equation with the activation energy 
Ai? ~ Ai^o — Fr, where Ai?o is the activation energy 
for r = and F is the linear response coefficient of the 
activation energy to the external torque t. 

In contrast to this, when the external torque is larger 
than a critical value, Z-form structure becomes unstable 
rather than metastable. For such a high external forces, 
the transition mechanism is switched to a spatially inho- 
mogeneous spinodal- like process (Sec. IIV C|) . Namely, the 
transition first occurs at the terminal which the external 
torque acts on. Then, the B-Z interface created there 
starts to sweep to the other terminal. When the inter- 
face reaches the other terminal, the whole transition is 
completed. In this case, the number of B-form segments 
scales as oc t^/^ with time t. We gave a phenomcnologi- 
cal estimate that accounts for this numerically obtained 
scaling, in which we assumed that the B-form DNA is 
approximated by a rigid rod and that the transition into 
B form, occurring at the Z-form segment adjacent to the 
B-Z interface, gives rise to the rotation of the B-form rod. 
Using these two assumptions, we reproduced the numer- 
ically obtained scaling relation as x/\/tI ~ [t/ naf'^'^ . 

We expect that the simple dynamical model developed 
in this paper equips the essence of interaction between 
DNA chiral structure and the mechanical response to the 
external forces, and hence we believe that the various 
structural transition phenomena elucidated in this paper 
can be verified in future experiments. 

It is noteworthy to remark that we have neglected 
the hydrodynamic interaction between DNA segments in 
modeling the effect of physiological bufi^er, as discussed 
in Sec. |lll This effect would become important especially 
for studying dynamical interaction between structures 
(e.g., B form, Z form) and three-dimensional configu- 
rations (e.g., plectoneme) of a DNA molecule, because 
two segments, far apart from each other in a nucleic acid 
sequence, can be spatially close to each other. The in- 
clusion of this effect is a further task of this work. 

As a final remark, we note that the dynamical aspects 
we have developed will provide a useful tool for investi- 
gating the dynamics of other structural transitions, such 
as a denaturation process [l^ US [HI . 



Appendix A: Kink propagation in a conserved 
system 

In this appendix, we show that kink propagations oc- 
cur in the conserved system (i.e, T = 0, k = 0). As 
an ideal case, we consider a DNA molecule of infinite 
length without external torques. From the principle of 
least action, we derive equations of motions for Q{s,t) 
and t). The action is given by 5 = Ldt, where L 
is the Lagrangian of Eq. The first variation of 5* is 
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given by 
dS 



{-h^ + c-^ix-e)}Sx]dsdt. (Al) 



Arbitrary Sx and SO should satisfy SS ^ 0, which gives 
the following equations of motion; 

(A2) 
(A3) 



/2^ + c^(x-e) = o 



ds 



To obtain kink solutions with velocities v, we set 

X(s, t) i^(s - ut), e(s, t) = G(s - vt) 

in Eqs. (|A2p and (|A3p . Thereby, these equations reduce 
to the following ordinary differential equations for F{x) 
and 0(0;) with x = s — vt: 



d^G d ,dV, 



da; 2 



dx"* dx dp' 



(A4) 
(A5) 



By eliminating F from these equations, the equation for 
G is given by 

d — A— - —( — ') -0 
^ dx'' da:^ dx dp 



with 



2 , -^2W2C 



(A6) 



This equation gives a first integral 



By setting g{x) 



cfG ^dG _dV_^ 
dx^ dx dp ^ 



^^^j— , one obtains 



J d^9 , dV^ . ^ 



This is the equation of motion for a particle of mass di 
moving in a potential of — V{g) — Cig. Hence, we 
obtain the following another first integral: 



According to Jensen et al. [l^, here we 

assume the following symmetric potential V(p) = 

r 9 i2 ^„ 
Vq {p/uj2) - 1) . For Eq. (gJI to allow for kink so- 
lutions, we have to choose Ci and C2 as follows: 



^0 f^o , A 



Vo V'^o 4 



The resulting equation is given by 



dg 
dX 



(A8) 



where a'^ ^ uj^ + ^ and X = xj^fdv^^JW^. This 
equation has kink solutions giX) = ±atanh(aX). 

Note here that the family of kink solutions exist for 
arbitrary v with < G/I2 [see Eq. (jA6P ]. which is in 
contrast to the result of nonconservative systems, where 
the speeds of interfaces were uniquely determined with 
the viscosity (Sec. lIVCp . 
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